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ABSTRACT 

We have modified the standard code for primordial nucleosynthesis to include 
the effect of finite nucleon mass on the weak-interaction rates as calculated 
by Seckel [p]]. We find a small, systematic increase in the 4 He yield, AY ~ 
0.0057 Y, which is insensitive to the value of the baryon-to-photon ratio r] 
and slightly larger than Seckel's estimate. The fractional changes in the 
abundances of D, 3 He, and 7 Li range from 0.08% to 3% for 10 -11 < rj < 10~ 8 . 



1 Introduction 



Primordial nucleosynthesis is one of the cornerstones of the hot big-bang 
cosmology. The agreement between its predictions for the abundances of 
D, 3 He, 4 He and 7 Li and their inferred primordial abundances provides its 
earliest, and perhaps most, stringent test. Further, big-bang nucleosynthesis 
has been used to provide the best determination of the baryon density || |f 
and to test particle-physics theories, e.g., the stringent limit to the number 
of light neutrino species [|] . 

The scrutiny of primordial nucleosynthesis, both on the theoretical side 
and on the observational side, has been constant: Reaction rates have been 
updated and the effect of their uncertainties quantified || , finite-temperature 
corrections have been taken into account |J, the effect of inhomogeneities in 
the baryon density explored 0, and the slight effect of the heating of neu- 
trinos by e ± annihilations has been computed ||; the primordial abundance 
of 7 Li has been put on a firm basis 0, the production and destruction of 
D and 3 He have been studied carefully |H|, and astrophysicists now argue 



about the third significant figure in the primordial He abundance |TT 



A measure of the progress in this endeavour is provided by the shrinking 
of the "concordance region" of parameter space. The predicted and measured 
primordial abundances agree provided: the baryon-to-photon ratio lies in the 
narrow interval 3 x 10~ 10 < rj < 4 x 10 -10 and the equivalent number of light 
neutrino species N v < 3.3 ||. The shrinking of the concordance interval 
motivates the study of smaller and smaller effects. 

Seckel QXJ] has recently calculated the corrections to the weak-interaction 
rates that arise from taking account of the finite nucleon mass (in the stan- 
dard code these rates are computed in the infinite- nucleon- mass limit). The 
corrections involve terms of order m e /m^, T/m^, and Q/m^, which are all 
of the order of 0.1%. Here m e is the electron mass, is the nucleon mass, 
T ~ 0( MeV) is the temperature during the epoch of nucleosynthesis, and 
Q = m n — m p = 1.293 MeV is the neutron-proton mass difference. The weak 
interaction rates govern the neutron-to-proton ratio and thereby are crucial 
to the outcome of nucleosynthesis; e.g., the mass fraction of 4 He produced 
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is roughly twice the neutron fraction at the time nucleosynthesis commences 
(T ~ 0.07 MeV). 

The net effect of the finite-nucleon-mass corrections is to decrease the 
weak rates by about 1% around the time of nucleosynthesis. Based upon a 
simple code that follows the neutron fraction Seckel estimated that the cumu- 
lative effect of all the corrections increase the mass fraction of 4 He synthesized 
by AF fa 0.0012. Because the third significant figure of the primordial 4 He 
abundance is now very relevant, we decided to incorporate the finite-nucleon- 
mass corrections to the weak-interaction rates into the standard nucleosyn- 



thesis code ||12|| . In the next Section we describe the modifications we made; 
we finish with a discussion of our results for the change in the yield of 4 He, 
which is slightly larger than Seckel's estimates, and for the changes in the 
yields of the other light elements. 



2 Modifications to the Standard Code 
2.1 Role of weak interactions 

The weak interactions that interconvert neutrons and protons, n <-> p + e + u, 
n + e <-> p + u, and n + v <-» p + e, play a crucial role as they govern the 
neutron fraction, and the neutron fraction ultimately determines the amount 
of nucleosynthesis that takes place. (Here and throughout we use e to indicate 
electron or positron, and v to indicate electron neutrino or antineutrino; 
the appropriate particle or antiparticle designation follows from charge and 
lepton number conservation.) 

The weak-interaction rate per nucleon is very roughly r n _> p ~ T p ^ n ~ 
G 2 F T 5 , while the expansion rate of the Universe H ~ T 2 /mp\, here Gp = 
1.1664 x 10~ 5 GeV -2 is the Fermi constant and m P1 = 1.22 x 10 19 GeV is the 
Planck mass. At temperatures greater than about 0.8 MeV, T is greater than 
H, and the neutron-to-proton ratio tracks its equilibrium value 

(-) = (m n /m P ) 3/2 e-<^. (1) 

\P / EQ 
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Two comments are in order. First, in the infinite- nucleon- mass limit the 
prefactor is unity; taking account of this factor tends to increase the equilib- 
rium ratio by about 0.2%, suggesting that the neutron abundance and final 
4 He abundance should be correspondingly larger. Second, at a temperature 
of about 1 MeV, the neutrino and photon temperatures begin to deviate as 
neutrinos decouple from the electromagnetic plasma (e ± , 7) and annihila- 
tions begin to heat the photons relative to the neutrinos as e ± pairs transfer 
their entropy to the photons. When this happens the equilibrium value of 
the neutron-to-proton ratio is no longer given by such a simple formula. 

When the temperature of the Universe drops below about 0.8 MeV the 
weak-interaction rate are no longer greater than the expansion rate and the 
neutron-to-proton ratio ceases to track its equilibrium value, and is said to 
"freeze out." Until nucleosynthesis begins in earnest (T ~ 0.07 MeV) the 
neutron-to-proton ratio decreases slowly due to weak interactions (especially 
neutron decay). For 77 ~ 3 x 10~ 10 it decreases from about 1/6 at freeze 
out to about 1/7 when nucleosynthesis begins, finally resulting in a mass 
fraction of 4 He equal to about 25%. For a detailed discussion of the physics 
of primordial nucleosynthesis see Refs. 



13, p. 



In the infinite-nucleon-mass limit the rates (per particle) for the six reac- 
tions that interconvert neutrons and protons are given by ]TB[ ] 

pE(E-Q) 2 dE 



p + e 



p + v 



X(n + v 
X(n + e 
X(n — * p + e + v 
X(p + e — > n + v 
X(p + v — > n + e 



X(p + e + v 
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, [ e (E-Q)/T v + !] [ e -B/T + !] ' 

P E(E + Q) 2 dE 
le [ e -(£+Q)/r„ + ij pTr + X ] 5 

Q P E(E -QfdE 

^ [ e &-Q)/T v + 1] [ e -E/T + !] ' 

pE{E -QfdE 

} [ e (Q- E )/ T » + 1] [e E ' T + 1] ' 

pE(E + Q) 2 dE 
le [e(s+Q)/r„ + 1 ][ e -B/r + 1 ]5 

Q pE(E-Q) 2 dE 



(2) 



[ e (Q-E)/T v + !] [gB/T + !] ' 

where T denotes the photon temperature, T u the neutrino temperature, and 
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the common factor 

p ^cos 2 ^(l + 34) 

C " 2^ /em - 

Here 9c is the Cabibbo angle (cos 0c = 0.975), ca = 1.257 is the ratio of 
the axial vector to vector coupling of the nucleon, and /em ~ 1.08 quantifies 
the radiative and Coulomb corrections, which for our purposes here are not 
important (for more details see Ref. 0). 



2.2 Finite- mass corrections 

Seckel has calculated the corrections to the weak rates due to finite nucleon 
mass 0. He has grouped them into three categories; the corrections due to: 
(i) thermal motion of the target nucleon; (ii) the recoil energy of the outgoing 
nucleon; and (iii) weak magnetism. The net effect of these corrections is to 
reduce the weak rates by about 1% around the time the neutron-to-proton 
ratio freezes out. This in turn causes freeze out to occur slightly earlier, at 
a higher value of the neutron-to-proton ratio, resulting in an increase in 4 He 
production. 

The corrections to the rates for neutron decay and inverse decay are the 
simplest, just the factor that accounts for time dilation, 

/ 3 T \ 

A(n <-> p + e + v) — ► ( 1 — ) X(n <-> p + e + v) . (3) 

Since the neutron lifetime is used to normalize all the weak rates, the recoil 
effect is automatically taken into account. 

The corrections to the weak interactions that involve 2 <-> 2 scatterings 
are organized in a way such that the integrands in the previous expressions, 
cf. Eqs. (H), are simply multiplied by a correction factor ,[] 

1 - <^n + 7wm + 7rec + 7th + 7etc, (4) 

where S n = —0.00201 "uncorrects" the neutron lifetime for recoil effects, 
which are taken into account in the term 7 rec , j wm is the weak magnetism 

1 Our notation differs very slightly from that of Ref. 7i = 7i/7o- 
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correction, 7th is the thermal correction, and 7 e t c is the sum of three smaller 
(by a factor of 100) and less important terms, which we have included in our 
computations, but which have negligible effect. 

The weak magnetism, recoil, and thermal corrections are given by |l[ 

2c A f 2 E^j + E^l 



7r 



1 + 3ci rriNEiEs 
7wm _ 2E 1 fc 2 + E 3 (£; 2 + fc 2 ) m^-m^-Q 2 

^ 2(1 + 34)77^3 

4 6E 2 E 3 - 6£iS| - 3Eikj - 4£ 3 £; 2 - E\k\ 



1 + 3c^ 2m N E 1 E 3 



7th 



/ T 



Vmjv 



3£ 2 + 2k\ 3kfE 1 + 3EfE 3 + 2£; 2 £ 3 k\E. 



2E 1 E 3 ' 2E x k\ 2k 4 
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(6) 
(7) 



subscript 1 (3) refers to the incoming (outgoing) lepton, subscript 2 (4) refers 
to the incoming (outgoing) nucleon, E{ is the energy of the zth particle, k{ 
is its momentum, and f 2 = ±3.62 where + applies to reactions with leptons 
and — to reactions with antileptons. For the three smaller terms embodied 
in 7 etc we refer the reader to Ref. 

Three small points; the limits of integration are those in the infinite- 
nucleon-mass limit, cf. Eqs. @, with one exception (see below). In Eqs. 
@ the integrals are performed over the energy of the electron; the neutrino 
energy is just that of the electron plus or minus Q. The integral of the 
final term in 7^, which is proportional to k% 4 , diverges for the reaction 
p+u — > n+e for the infinite-nucleon-mass limits because k 3 —>■ at threshold. 
For this term, and only this term, the lower limit of integration involves 
the finite-nucleon-mass kinematics: the minimum electron momentum is: 
(&3)min = (^e/^v) 1//2 (Q + 7™ e ) [0 • Because essentially all of the integral 



accumulates near (k 3 ) m i n this term can be integrated analytically: 

P E(E + Q) 2 dE k\El _ (T/m N )(Q + m e fm z J 2 m^ 2 



m N J J [ e ( E +Q)/ T » + l}[e~ E / T + 1] 2fcf 2[e^ +m ^/ T " + l}[e~ m ^ T + 1] ' 

(8) 

Finally, we mention that Seckel H has derived linear fits to the perturbed 
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weak-interaction scattering rates: 



A(n — > p) 
\{p — > n) 



1 -<L -0.00185 - 0.01032 



1 - 5„ + 0.00136 - 0.01067 



T 

m N 

(- 



\{n^p); (9) 
A(p^n); (10) 



where the first linear correction applies to the 2 2 reactions that convert 
neutrons to protons, n + v — > p + e and n + e — > p + z/, and is valid for 
2MeV > T > 0.3 MeV, and the second applies to the 2 <-> 2 reactions that 
convert protons to neutrons, p + v — > n + e and p + e — > n + u, and is valid 
for 2 MeV > T > 0.7 MeV. 



3 Results and Conclusions 

We have modified the integrands in the subrountines that compute the weak- 



interaction rates in the standard nucleosynthesis code |L2] as outlined above. 
For comparison, we have also modified the standard code just using Seckel's 
linear fit, which is much easier to implement since the usual rate is multiplied 
by a factor outside the integral. Our results, which were obtained by taking 
three massless neutrino species and a mean neutron lifetime of 889 sec, are 
shown in Figs. 1 and 2. 

In Fig. 1 the change in the mass fraction of 4 He and 0.0057 times the 
4 He yield are shown. While over the range rj = 10~ n to 10 -8 AY varies 
from about 0.0004 to 0.0015, AY is remarkably close to 0.0057 V. We also 
computed the change in 4 He yield using Seckel's linear fit; for the values of 
rj above, AY/Y was slightly lower, by about 0.05%. Neglecting the subdom- 
inant terms (7 e tc) decreases AY/Y from about 0.57% to about 0.53%. 

In Fig. 2 we show the fractional changes in the abundances of D, 3 He, 
and 7 Li (relative to H). These changes, over the same interval in baryon-to- 
photon ratio, range from 0.08% to almost 3%. Since the inferred primeval 
abundances of these elements are no where near as well known as that of 
4 He, these changes are of little relevance at present. 
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FIGURE CAPTIONS 



Figure 1: The change in the yield of 4 He as a function of the baryon-to- 
photon ratio r], and, for comparison, 0.0057 times the 4 He mass fraction. 

Figure 2: The fractional changes in the yields of D (solid curve), 3 He (broken 
curve), and 7 Li (dotted curve) as a function of the baryon-to-photon ratio r\. 
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